Mecânica estatística 


quântica II: bósons 


META 
Construir o formalismo da mecânica 


estatística quântica para bósons. 


OBJETIVOS 

Ao fim da aula os alunos deverão ser 
capazes de: 

Entender a influência das propriedades 
quânticas na construção do formalismo da 
mecânica estatística. 

Obter a distribuição de Bose-Einstein. 
Abordar a radiação do corpo negro e o 
calor específico de sólidos. 

Resolver problemas envolvendo estes 


conceitos. 


PRÉ-REQUISITOS 
Aulas anteriores e cálculo diferencial e 


integral. 


*Bósons são partícu- 
las quânticas que 
possuem número quân- 
tico de spin inteiro 
(s=1,2,3,...) como, 


por exemplo, o fóton 
(s=1). Estas partícu- 
las possuem função de 
onda simétrica e não 
satisfazem o princípio 
de exclusão de Pauli 


[1]. 


(160, 
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9.1 Introdução 


Caro aluno, como foi visto na aula anterior, ao levarmos em 
conta as propriedades quânticas das partículas, estabelecemos re- 
gras de contagem que diferem do caso clássico, dando forma a uma 
estatística distinta da de Maxwell-Boltzmann. Vimos que ao levar 
em conta o princípio de exclusão de Pauli (PEP) para férmions, 
obtemos a distribuição de Fermi-Dirac. No entanto, como tratar 
os bósons estatisticamente? Estes são partículas quânticas que não 
respeitam o PEP, pois a função de onda de bósons é simétrica com 
relação a troca de partículas. No entanto, assim como os férmions, 
bóson são partículas indistinguíveis por consequencia do princí- 
pio da incerteza de Heisenberg, que é uma das maiores bases da 
mecânica quântica. 

Nesta aula, iremos reformular a mecânica estatística com base 


nas propriedades quânticas dos bósons. 


9.2 Distribuição de Bose-Einstein 


Apesar dos bósons não respeitarem o PEP, estes devem satis- 
fazer um princípio mais fundamental ainda: a indistinguibilidade. 
Com base nesta característica, vamos encontrar a distribuição que 


rege os sistemas bosônicos. 


9.2.1 Contando microestados 


Considere, por simplicidade, um sistema com N bósons não in- 
teragentes e distribuídos num espetro de energia e € (gH- Cada 
valor de energia e; possui n; bósons em 9; níveis, ou seja, estados 


com degenerescência g;. Neste caso, a contagem de microestados 


Introdução à Física Estatística 
ececjo/00 ||00 | 0 ce) 0 jo 0) 6 


Figura 9.1: Uma possível arrumação de n; = 16 bósons (bolas escuras) 
em 9; = 9 níveis, os quais estão esquematicamente separados na figura por 


9; — 1 divisórias (linhas verticais). (Figura adaptada da ref. [1].) 


é equivalente à contar o número de maneiras em que n; bolinhas 
indistinguíveis podem ser distribuídas em 9; caixinhas. Esquemati- 
camente, isto é equivalente a colocar as n; bolinhas numa fileira e 
separá-las com 9; — 1 divisórias, como mostra a fig. 9.1 , onde ve- 
mos um caso particular de n; = 16 bolinhas e g; = 9 caixinhas. Da 
esquerda para direita, vemos três bolinhas na primeira caixinha, 
uma bolinha na segunda caixinha, etc. Sendo assim, o número 
de maneiras distinguíveis que podemos distribuir as bolinhas nas 
caixas é igual ao número de permutações possíveis de n; + 9g;— 1 
elementos, sendo n; bolinhas e g; — 1 divisórias. Este número de 
permutações corresponde a (n; + 9; — 1)!. No entanto, a permu- 
tação do número de bolinhas (n;! possibilidades) e do número de 
divisórias [(9; — 1)! possibilidades] não altera o microestado físico. 


Portanto, concluímos que há 


(n; + 9; — 1)! 
n!(g; — 1)! 


microestados para a energia e;. Sendo assim, obtemos o número 


total de microestados 


L 
w(m,..,n)= | (utah (9.1) 
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9.2.2 Maximizando o número de microestados 


Seguindo a prescrição da sec. 5.4.2 , fazemos a maximização de 


W e obtemos 


n; = E (9.2) 


onde continuamos com 
1 


Pes (9.3) 


A probabilidade de um bóson ter energia e;, é obtida fazendo 


P;=nj/N leg. (5.7)| e, sendo assim, temos 


Pp; = “ : 
7 Nexp(a+e;/kBT)-—1 


(9.4) 

A constante a é encontrada com a condição de normalização 
da probabilidade, a P; = 1. Note que, como os bósons não 
seguem o PEP, não faz sentido definirmos a constante a em ter- 
mos da energia de Fermi. No entanto, como a probabilidade deve 
sempre ser positiva, entre O e 1, a constante a deve ser não nega- 
tiva, ou seja, 


a>o0. (9.5) 


9.3 Radiação do corpo negro: gás de fótons 


O fóton é um bóson, que corresponde à menor divisão de uma 
onda eletromagnética. Pela teoria da relatividade restrita, a ener- 


gia de uma partícula é [1] 


Ei ae fail DI 2 
e=me = a/méc +p”, 


onde c é a velocidade da luz e 


mo =mvy'1— v2/c? 
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é a massa de repouso da partícula. No entanto, para o fóton mo = 


0. Sendo assim, a energia do fóton é 
€= cp. 


Em mecânica quântica, podemos usar a definição de operador 


diferencial do momento [2], 


d 
EN ad 
p= 


para resolver a equação de Schródinger de um fóton em uma caixa 


cúbica de aresta a e, com isso, obter 
(9.6) 


onde k = VR+R+ e 31,J2,)3 € N. Podemos considerar que 
a caixa é muito grande, ou seja, a >> À, onde À é o comprimento 
de onda do fóton. Sendo assim, como o espaçamento de níveis de 
energia é da ordem de ch/a [eg. (9.6)|, podemos tomar o limite do 
contínuo e obter a densidade de níveis de maneira análoga ao que 


fizemos na sec. 8.3* 


g(e) = 23 fa kó (e — et) ; (9.7) 
Resolvendo esta integral, obtemos 
81V 
g(e) = TE é. (9.8) 


No limite do contínuo, podemos reescrever a eq. (9.4) da seguinte 


forma 
= dl g(e) 
RlejES Nexp(a+e/kBT) — 1 


Agora, podemos considerar que a caixa cúbica que confina os fó- 


tons em nosso problema é uma cavidade que absorve e reemite 


“O fator multiplicativo 
2 visto na sec. 8.3 
por consequência 
do spin do elétron, 
continua sendo válido 


para fótons, pois 
estes possuem duas 
direções de polarização 
independentes (ver- 
tical ou horizontal), 
visto que as ondas 
eletromagnéticas são 
transversais. 


(163, 


9— B 


(164, 
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fótons até alcançar o equilíbrio, quando os números de fótons ab- 
sorvidos e reemitidos são iguais. Este é justamente o problema da 


radiação do corpo negro. Como neste caso o número de fótons não 


é constante, temos dn; + 0, implicando em a = O (ver ativ. 


9.1). Logo, 


e g(e) 
a Nexp(e/kBT)—1 
8tV é? 


Nc'h3 exp(e/kBT) — 1º (9:9) 


Quanticamente, podemos escrever a energia do fóton da seguinte 
forma 


e=hy, 


onde v é a frequência do fóton. Sendo assim, é conveniente calcular 


a distribuição de frequências da seguinte forma 
ú(v)dv = p(ejde = p(hv)hdv. 


Logo, concluímos que 


8tV e: 
NS co (85) 
Aqui, N pode ser calculado através da normalização de densidade 


de probabilidade 
d Ú(v)dv = 1. 
0 


Também podemos calcular a densidade de energia por volume 


Ú(v) = 


(9.10) 


em função da frequência, fazendo 


hv 
Elm = Ni 90). 


Sendo assim, obtemos 


e(v)= (9.11) 
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conhecida como lei de radiação de Planck, a qual foi sugerida por 


Planck para o problema da radiação do corpo negro [1]. 


Exemplo 9.3.1. Discuta as transições atômicas espontâneas e in- 
duzidas por radiação num sistema de dois níveis. (Exemplo adap- 


tado da ref. [1].) 


Solução: Tomamos como base a fig. 9.2 onde temosN; (No) áto- 
mos com energia Fi(E5) e Es > E, com E — E, = hv. 
Os átomos que estão no nível E» podem passar espontanea- 
mente para o nível E1, com probabilidade por unidade de 
tempo As1. Considere que há radiação com frequência v e 
densidade de energia €(v), interagindo com os átomos nos 
dois níveis de energia induzindo transições por absorção ou 
por emissão. É razoável supor que a probabilidade transição 
induzida por radiação é proporcional a E(v). Sendo assim, 
podemos dizer que a probabilidade por unidade de tempo 
de transição induzida por radiação é B»€(v), do nível E» 
para o F1, e B12€(v), do nível E; para o E5. Assim, con- 
siderando a absorção e emissão de radiação dos átomos no 
nível Es, podemos estimar a variação de átomos neste nível 


por tempo da seguinte forma: 


dN: 
E = Bé(v)Ni — [An + Bné(v)]No. 
e E O rg A 


absorção emissão 


Quando o equilíbrio equilíbrio entre os átomos e a radiação 


é estabelecido, devemos ter dNo/dt = 0, ou seja, 


BE(v)N = An + BaE(v)No. 


(166, 
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N5 átomos 
E, 


Am + Bn€(v) l Í B,2€(v) 
É1 


N, átomos 


Figura 9.2: Transições espontâneas e induzidas. (Figura adaptada da ref. 


[1]. 


Se os átomos estão em equilíbrio térmico e seguem a estatís- 


tica de Maxwell-Boltzmann, temos 


Logo, 
hv 
B,€(v) exp | — | — An + Bne(v) 
kBT 
e, consequentemente, 


An/Bio 


Sus 
exp (Er) — Ba/Bio 


Comparando esta equação com a eg. (9.11), obtemos 


p= pet 1 
Bio ER Bio 
Com isso, podemos ver que 
prob. de emissão espontânea Am o hv o 
prob. de emissão induzida Ba€(v) - kBT , 


Perceba que para hv >> kpT' a emissão espontânea é muito 
maior que a induzida, a qual pode ser desprezada. No en- 
tanto, para hv < kpT, a probabilidade de emissão induzida 


é muito maior que a espontânea. 
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9.4 Calor específico de sólidos: gás de fônons 


As vibrações coletivas de átomos em sólidos são responsáveis 
pela variação de energia do sistema através da variação da tem- 
peratura. Sendo assim, iremos considerar estes efeitos vibracionais 
para calcular o calor específico dos sólidos. 

Estas vibrações coletivas formam ondas estacionárias, que po- 
dem ser interpretadas como ondas elásticas. As frequências destas 
ondas estacionárias dependem da forma e do tamanho do sólido. 
Existe um bóson para estas ondas equivalente ao fóton para a onda 
eletromagnética, o qual é chamado de fônon. Este nome é justifi- 
cado pela natureza destas ondas elásticas, já que se propagam com 
a velocidade do som. Podemos, portanto, modelar a vibração da 
rede cristalina de um sólido através de gás de fônons, de maneira 
análoga ao que fizemos na seção anterior, onde modelamos o prob- 
lema da radiação do corpo negro como uma gás de fótons. 

As ondas elásticas em um sólido possuem duas classes, longitu- 
dinais e transversais, as quais se propagam com, respectivamente, 
velocidades v; e v;. No entanto, as ondas longitudinais possuem 
apenas um único grau de liberdade, enquanto as transversais pos- 
suem dois. Sendo assim, podemos manter a densidade de níveis da 
eg. (9.8), trocando c por v; ou v, e dividindo por 2 para as ondas 


longitudinais. Logo, 


81V 5 
gi(e) = EE ; 
47V » 
gi(e) = ESTE : 


Sendo assim, a densidade total de níveis é 


o 2 4V/1 Ag 
s(e) = gil) + a(o) = E (5 + 5) E. 
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Em um meio contínuo não há limite para o número total de modos 
de vibração, ou seja, o número de níveis. No entanto, em um 
sólido com N átomos, todo modo vibracional é descrito por 3N 
coordenadas espaciais dos átomos, sendo este o número total de 
modos independentes de vibração. Com isso, devemos impor o 
limite superior de energia, hvo, o qual define a frequência de corte 


vo. Com isso, temos 


a ARM fd DO 
3N - [ g(ejde = ER (a + 5) À ede, 


e, portanto, obtemos uma relação que determina a frequência de 


9N 1 2 
vg V; Vs 


corte 


Logo, podemos escrever a densidade de níveis da seguinte forma 


9N » 
g(e) — n3vÊ E 


Com isso, podemos escrever a distribuição de energia do fônon 


dO) = rã 
M exp(e/kBT) — 1 
9N é 
Mh3vi exp(e/kBT) — 1º 


(9.12) 


onde M é o número de fônons em equilíbrio térmico com o sólido. 
Análogo ao fóton, podemos escrever a energia do fônon como 
e = hv onde v é a sua frequência. Agora, podemos calcular a 


distribuição de frequências da seguinte forma 


onde obtemos 


(9.13) 
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C (graphite) 
O AlOg 
x KCl 


0 0.5 1.0 1.5 2.0 FIN 


Figura 9.3: Calor específico de uma rede cristalina em função da temper- 


atura. (Figura retirada da ref. [1].) 


Note que, ao impor a normalização desta densidade de probabili- 

dade, é possível obter o número de fônons, M, o qual dependerá 

da temperatura, do número de átomos e da frequência de corte. 
Agora podemos calcular a energia interna vibracional do sólido 


fazendo* 
vo 
U=M(e=Mh(v) = Mi [ vo(v)dv, 
0 


implicando em 


(9.14) 


Usando a eq. (7.1), podemos obter o calor específico molar a 
volume constante do sólido 


3 OpD/T 4 
PRA e (55) / RE qdo 
0 


ni (9.15) 


onde 


“Poderíamos adicionar 
a energia de ponto 
zero associada a cada 
modo de vibração, mas 
esta não interfere na 
obtenção do calor es- 
pecífico, pois não de- 
pende de temperatura. 
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é a temperatura característica de Debye. Como vo depende do ma- 
terial, a temperatura de Debye também é característica específica 
deste. 

Resolvendo numericamente esta integral, alcança-se o resultado 
apresentado na fig. 9.3 , o qual concorda muito bem para vários 
materiais. O regime assintótico de cy em altas temperaturas é 
constante e igual a 3R. Este comportamento é conhecido como a 


lei de Dulong-Petit. 


9.5 Conclusão 


Vimos que ao considerarmos as propriedades quânticas de in- 
distinguibilidade para os bósons, obtemos a distribuição de Bose- 
Einstein, a qual é fundamental para estudarmos estatisticamente 
sistemas bosônicos. Como exemplo, abordamos o gás de bósons 
modelando a radiação do corpo negro e o calor específico de sóli- 


dos. 


9.6 Resumo 


A distribuição de Bose-Einstein para bósons é dada por 


Pp; = É 
7 Nexp(a+e;/kBT)-—1 


Para o problema da radiação do corpo negro, modelado por um 


gás de fótons, a distribuição de frequência para os fótons é 


O número de fótons em equilíbrio com as paredes da cavidade, 


N, pode ser calculado através da normalização de densidade de 
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probabilidade 
/ úÚ(v)dv = 1. 
0 


O calor específico em um sólido é ocasionado basicamente pela 
vibração da sua rede cristalina. O fônon é o bóson associado a 
esta vibração. Modelando este problema como um gás de fônon, 


obtemos a distribuição de frequência dos fônons 


— 9N v2 
— Mi h 
V0 exp (dr) —1 


Ú(v) 


Note que, ao impor a normalização desta densidade de probabili- 
dade, é possível obter o número de fônons, M, o qual dependerá da 
temperatura, do número de átomos e da frequência de corte. Com 


este resultado, obtemos o calor específico molar a volume constante 


do sólido 
TN3 pOD/T 44 
cy =9R (55) i asno) 342, 
Op/ Jo — lexp(z)-1] 
onde 
hvo 
op="0 
D RE 


é a temperatura característica de Debye. 


9.7 Atividades 


ATIV. 9.1. (a) Mostre que ao maximizar o número de microes- 
tados da eq. (9.1) se obtém a eg. (9.2). (b) Usando os passos da 
resolução do item anterior, justifique o motivo pelo qual a = 0 se 


o número de bósons não se conserva. 
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Comentário: Esta maximização pode ser feita com o uso de mul- 
tiplicadores de Lagrange. Uma consulta na ref. [1] poderá 


auxiliar na resolução desta atividade. 


ATIV. 9.2. Calcule a integral da eq. (9.7) até obter a densidade 


de níveis da eq. (9.8) para o gás ideal de fótons. 


Comentário: Esta integral pode ser resolvida como feito na sec. 
5.5. Lembre-se que o volume de um cubo é a sua aresta ao 


cubo. 


ATIV. 9.3. Considerando o problema da radiação do corpo negro, 
(a) mostre que o número de fótons dentro da cavidade em equilíbrio 


com os átomos das suas paredes é 


167V 
ds h3c3 


((3)(kBT)º, 


onde a 
do = DE 
j=1 


é a função Zeta de Riemann. (b) Calcule a densidade de fótons na 
cavidade a 300 K. 
Dados: (3) = 1,202056903. 


Comentário: Considere a normalização da densidade de proba- 
bilidade da eq. (9.10). A resolução da integral pode ser sim- 
plificada se multiplicarmos o numerador e o denominador do 
integrando por exp(—hv/kpBT). Em seguida identifique o fa- 


tor deste integrando que corresponde à soma de uma P.G. 
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infinita e convergente. Com isso, é possível transformar a in- 
tegral em uma soma de infinitas integrais, as quais podem ser 
facilmente resolvidas com técnicas de derivação paramétrica 


que você já está habituado. 


ATIV. 9.4. Para o problema da radiação do corpo negro, (a) 


mostre que o n-ésimo momento da frequência é 


mo (N+2BCn+3) E 
tis o) VEDA 


(b) Com este resultado, mostre que a energia média (energia in- 


terna) por fóton e sua variância são 


4 
(e) = 30 7 er 27kB7, 


6(8) 


S(3) 
(c) Interprete estatisticamente o resultado do item anterior. 


Dados: €(3) = 1,202056903; €(4) = 14/90; €(5) = 1,036927755. 


de= vvar(e) = [es — 9 El “ksT a 1,75kBT 


Comentário: Lembre-se da definição de momentos na estatística 
abordado na sec. 4.4 . Recomenda-se fazer a atividade ante- 


rior antes desta. 


ATIV. 9.5. Estime a razão entre as probabilidades de emissão 
espontânea e induzida a 300 K para (a) a região de micro-ondas, 
v = 10% Hz, e para (b) a região ótica, v — 105 Hz. (Atividade 


adaptada da ref. [1].) 
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Comentário: Use o resultado do exemplo 9.3.1. 


ATIV. 9.6. Com base no exemplo 9.3.1 , justifique o fato de que 
as radiações emitidas de forma induzida são coerentes e de forma 


espontânea são incoerentes. 


Comentário: Uma consulta na ref. [1] é muito recomendada para 


a resolução desta atividade. 


ATIV. 9.7. (a) Demonstre a eq. (9.15) partindo da eg. (9.14). (b) 
Com o resultado do item anterior mostre que cy = 3R para altas 
temperaturas (7 >> Op), o que corresponde à lei de Dulong-Petit. 
(c) Justifique o resultado do item anterior com base no princípio de 
equipartição da energia. (d) Mostre que em baixas temperaturas 


(T < Op), 0 calor específico segue a forma 


(e) Com base no resultado do item anterior, determine (quantita- 
tivamente) uma maneira de se estimar experimentalmente a tem- 


peratura de Debye para um sólido. 


Comentário: Uma consulta na ref. [1| poderá auxiliar na res- 
olução dos itens (a) e (c) desta atividade. No item (b), como 
a temperatura é alta, é necessário aproximar o integrando 
da eg. (9.15) para pequenos valores de x e só em seguida 
fazer a integração. Note que, no item (d), como a temper- 


atura é baixa, é possível tomar o limite superior da integral 
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da eg. (9.15) como infinito, resultando numa integral muito 
semelhante às resolvidas nas atividades 9.3 e 9.4. Um esboço 
de pontos experimentais de cy em função de Tº para baixas 


temperaturas poderá auxiliar a resolução do item (e). 


ATIV. 9.8. Discuta a relevância do calor específico de um sólido 
provocado pela movimentação dos elétrons comparado com a con- 


tribuição da vibração dos átomos da rede cristalina. 


Comentário: Uma consulta na ref. [1] é muito recomendada para 


a resolução desta atividade. 


9.8 Próxima aula 


Na próxima aula, abordaremos o gás ideal quântico e com- 
pararemos as três estatísticas apresentadas até aqui: Maxwell- 


Boltzmann, Fermi-Dirac e Bose-Einstein. 
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